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maximum  rigid  body  acceleration  of  the  projectile 
material  constants,  eq.  (8) 
constants  due  to  body  forces 
dilatation 

normal  strain  components 

Young’s  modulus  of  elasticity 

body  force  In  the  z  direction 

coefficients  of  eqs.  (21)  from  Raylelgh-Ritz  method 

shear  modulus 

length  of  propellant 

radius.  Inner  radius,  outer  radius 

total  potential  energy 

displacement  along  r  axis 

elastic  strain  energy 

potential  energy  of  external  forces 

displacement  along  z  axis 

Rayleigh  Rltz  undetermined  coefficients 

shear  strain  component 

weight  per  unit  volume 

normal  stress  components 

shear  stress  component 
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INTRODUCTION 


The  launch  phase  of  a  gun-launched  rocket  is  associated  with  high 
acceleration  loading  of  up  to  10,000  g.  As  a  result  of  the  high  accel¬ 
eration  the  motor  case  portion  of  the  rocket,  which  contains  the 
propellant,  is  subjected  to  axial  and  lateral  loads.  The  lateral  loads 
are  associated  with  the  deformation  of  the  contained  propellant.  This 
"slumping  propellant"  behavior  is  the  concern  of  the  present  investi¬ 


gation. 


ANALYSIS 


The  propellant-motor  case  system  may  be  Idealized  as  a  two-dimensional 
system  because  of  axial  symmetry,  with  cylindrical  coordinates,  r  and  z, 
a  typical  section  is  shown  below. 


r  o - > 


fore  /^propellant 

r  *  a  - r — 


n 


■motor  case 


2.  w 


aft 


r,  u 


(ro»z)  the  resultant  boundary  conditions. 


u  *  aispiacement  in  the  r  direction 
w  »  displacement  In  the  2  direction 
r1*  Inside  radius  of  propellant 
rQ=  outside  radius  of  propellant 
L  *  length  of  propellant 
The  particular  problem  of  the  unbonded 
propellant  Is  considered  here.  The 
stiff  motor  case  shell  is  assumed  to 
completely  restrain  axial  displacement 
at  the  aft  end  (r ,0)  as  well  as  radial 
displacement  along  the  outer  radius 


w(r,  0) 
u(r0.z) 


0 

0 


crz 


(r ,0)  > 

(r0.z) 


lrz'ro 


The  resuming  sides  (r,  L)  and  (r,,z)  are  stress  free,  that  Is 
»z(r.L)  ■  0  Trz(r,L)  -  0 
or<ri.T>*0  Vz(rcz)  -0 


(la) 

(lb) 

(lc) 

(ld) 


The  acceleration  of  the  propellant  subjects  every  particle  to  the  body 
force 

fz  *  "pam  (2) 

where  p  is  the  weight/unit  volume,  and  am  is  the  maximum  acceleration 
In  g's. 
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A  solution  to  the  problem  is  obtained  by  the  Raylelgh-Rltz  method. 
According  to  this  method,  an  approximate  solution  satisfying  the  displace¬ 
ment  boundary  conditions  is  formed.  The  unknown  coefficients  of  the  solution 
are  determined  by  minimizing  the  total  potential  energy  of  the  system. 

Here  the  displacement  fields 

u(r,z)  =  (L-z ) (a  +  sr  +  yr2  +  6r3)  =  (L-z)f(r)  (3) 

w(r,z)  *  z(L  -  j)(t\  +  xr  +  cr2  +  ?r3)  =  z(L  -  |-)g(r) 

are  assumed.  The  linear  z  term  In  the  u  displacement  and  the  quadratic  z 
term  in  the  w  displacement  were  chosen  to  satisfy  the  following  conditions: 
i)  Results  from  the  Rohm  and  Haas  finite  element  solution, 
ref.  1,  Appendix  B,  for  the  accelerated,  unbonded  propellant 
show  or  and  oz  to  be  linear  functions  of  z.  Note  that  Equations  (3) 
give  this  condition,  i.e. 


11)  The  propellant  Is  assumed  to  be  nearly  Incompressible,  and  thus 
linearly  varying  compression  along  the  z  coordinate  causes  the  pro¬ 
pellant  to  move  towards  the  centerline  since  It  Is  restrained  by  the 
case  from  moving  outward,  l.e.  u  <;  0  for  all  z,and  u  =  0  for  z  *  L. 

This  assumed  displacement  field  Is  not  valid  for  the  bonded  propellant. 

The  equation  for  w  satisfies  the  displacement  boundary  condition  w(r,0)  *  0. 

The  displacement  boundary  condition,  u(rQ,z)  =  0  is  satisfied  by  setting 

a  «  -Br0  -  yr*  -  5ro  (4) 

The  total  potential  energy  of  the  system,  T,  resulting  from  the  displace¬ 
ment  fields  given  by  equations  (3)  must  be  formed.  Denoting  the  elastic  strain 
energy  as  U  and  the  potential  energy  of  external  forces  as  V,  we  have 

T  «  U  +  V  (5) 
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Determination  of  the  elastic  strain  energy  U  proceeds  as  follows. 


According  to  the  strain-displacement  relations 
ef  *  |p  ■  (L  -  z) (0  +  2yr  +  36r2) 

e6  =  £  =  (L  -  z)(£+  5  t  Yr  +  6r2) 

ez  *  fz“  c  (L  -  z)(n  +  xr  +  zr*  +  ?r3) 

Yrz  “  fz  +  fr  a  ”^a  +  0r  +  Yf2  +  6r^ 

!  2 

+z(^  -  z)(x  +  2cr  +  3cr  ) 

and  the  stress-strain  relations, 

or  s  Aef  +  B(ez  +  e0) 

o  *  Ae„  +  B(e  +  e  ) 

r  e  '  r  z1 

oz  =  Aez  +  B(er  +  efi) 

t  *  Gy 

rz  ’rz 

where 

B  ’  |T  V  vT(T  -  7vJ 

and 

6“  7(r+  vr 


(6) 


(7) 


(8) 


(9) 
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Substituting  equations  (6)  Into  (7)  yields  stress-displacement 
equations.  The  elastic  strain  energy  of  the  system  Is  given  by 


°1j  «1j  dv  "  I  III  (°rer  *  Ve  +  °zez  +  VzVz'  dv  (,0) 


u.I 


where  the  Integration  Is  over  the  volume  of  the  propellant.  The 
Independence  with  the  e  coordinate  gives 

r1  L 

u  *  2’  '  7  ||  (arer  +  aees  +  °zez  +  TrzVz>  r  dr  dz 

ro° 

For  convenience,  let 


(11) 


||  °rer  r  dr  dz 


II 


IT  1 1  a0ee  r  dr  dz 


(12) 


then 


||  °zez  r  dr  dz 
||  TrzYrz  r  dr  dz 


u  *  +  u.  +  u„  +  u,. 

abed 


(13) 


Substituting  the  stress-displacement,  and  strain  displacement  relations 
Into  equation  (12),  and  Integrating,  we  obtain 
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-a  -  T 


*h 

+  B  {^/i  +  £  (« n  tyi1  nr*)  + 1  (  ,jjj> 

+Zt\ 4-3<JoC  -hySX  )  4-  (4-yflJ  >A?fX+34^ 

+  a.n+^0  *4?  (syJ  +  aSX+ySq  rt 

+  3.-^)  +i‘  (3i1+si<+a»c)  +  li?A7|  * 


Jb  *  T 


|a  /i  *  £  ( jsN-a**^ 

+  |.7r5^8  +  -LjVJ 

-?-B[A(o<Tl4-o(j8)  +  ^  U^jr+ya^-i-yaVotX) 

+  *?{2<S+*flf  +  Y*l  +  0X4-c <5) 

+  4.,(<VySi  +  ayltSn+TA  a-(8f  +  .<«;) 

+  i’(sri  v^5)+a.‘(Si>+i,+ri;) 

*“C 


_  wu 

Jc  a"y 


[A{-£  Yl^+T  ^  aA  +  +a>,l*)+  A?  (a.  As  +3l>|4) 

4l!’(fVaX«)  +  i  $c*’+ 

+  B  j^V1  +£‘(4^>t  +  tX)  4  y  +*$) 

+£  (4  4*l 4  m  +  €  ■*■  «*  c)+-  4s  (  4-  JA  -v  32Pf 

+  a^)  +l‘(+4?43rc)4  ii^A7}!*’ 
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Ud  S  "G  jj-j^o(V*+f  <*/}»?  +  4  A«r)+  yUocS 

+  A*yfl)  +  ^£$/lfc  +  J  3<JA7} 

“  ucH00^* + P*)  +^(3*5  +*p% 

+^x)  +  j£  (3y8C  +  a*'§  4-iA)+  lfc(3tc 
+  Afi)  +}i?A7} 

+rolixV  +  fx§A3+A4(^+^^)+T^Al  (17) 

The  potential  energy  of  external  forces  V,  associated  with  the  Inertia 
body  force  f  given  by  equation  (2)  Is, 

ro  L 

V  *  -  |  fz  *  w  dV  *  -  2ir  p  am  *  w  *  r  dr  dz  (18) 

r1  0 

Substituting  the  expression  for  w  from  equation  (3)  Into  equation 
(18)  yields 

3  r 

_|  2931495° 

V  *  -  -j-  p  am  '  [nr  +  j  xr  +  j  Cr  +  g  tr  ]  (19) 

r1 

The  total  potential  energy  of  the  system  is  obtained  as 

T  ■  vvw»  <20> 
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The  thee, era  of  minimum  potential  energy  states  that  equilibrium  Is 
associated  with  a  minimum  of  T.  Here  T  Is  a  function  of  the 
Independent  parameters  B.y.d.n.e  and  The  minimum  of  T  Is  obtained  by 
taking  the  partial  derivative  of  T  with  respect  to  each  of  these  seven 

coefficients.  This  yields  the  system  of  7  linear  algebraic  equations 
In  7  unknowns, 


F1Z  ^12  •  • 

•  •  •  •  Fj  7 

6 

c; 

^21  ^22  •  • 

•  •  •  •  F27 

Y 

c2 

• 

• 

6 

• 

• 

• 

n 

’ 

• 

• 

• 

X 

• 

• 

• 

C 

• 

F71  F72 

F77 

c 

4 

C7 

The  explicit  expressions  for  the  F(j  and  C,  (1(J  .  1,  7)  coefficients 
are 


F(',0  =  A  {-a*;-™?  +4*. A;  +  JU  Af,  } 

+B  {-ahj-aa?}  +  3G  { n? aJ  1 4  „.*>  - 1  a*  j  LA 

=  M-^*-***  4A.at+JA(,JU  A,  j 

+  +  +  3G/V2  »*-«*,* 

^  nxni  .  4.  _V  Ai  130  * 

-t|r  Ai  4-  4*  -  |  4*  l 


F(a,»)  «  F(i,*) 
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fM~f(3j)=  ii.Aj  +  aa,-U£} 

+  vaA./\J  -a*? 

+  3<*£ [1  ni-  Jli  A?  +  fnl*\  +  jfloA*  -  j-*‘  } 

F(i,4)*F(4j)  =  B{-AAi +aA0Ai\ 

F(iJs)  =  F(%.)=  e'.i(/lo-Ai)-JlJ  +  ».<x  \ 

-  JG^-±-  *2  +  i  *0*1-  s*l\ 


fO>}«fM= 


=  b{|aJ  +  ja,a\-a1  } 

'  aM"t  /i4  + 1  n,*l  -  i  a*  } 


F(|,t)-F(74i)-  8{i*:+i*.A!-  |,n 

-*1-*^ +**•<- hi 

F(a,i)  -  A  {‘t-*.4-  i**  +  * J»i  Aj-  +  a/I,4  in  J2 1 
4B|-2aJ  +  4<  /£-**?} 

+  3&/^{  -  aJa'  +  aJA;-^  A;  | 

F  (a,s)  *  F  (3,*) «  A  {  g  a?  -  £  iij  + 1  a  *  +  *i  A;  af  U.  As  ] 

+  B  {-aaJ  +a*{  a?  +aa^  a*~  a  a*} 

♦»  ig  «I-  vs  «M<:  a?  +  |4a*-  i  < } 

F(a,4)  =  F(4,a)  =  +  1/l^Ai} 

FM*F(S,x)  =  B\i(/l4-Ji4)-/t'+A,‘Ai| 

-a&|- ^  A,  +  •£  <2  A*  - 4*  } 


_ — 


P(a,0-Fl%a)*  Blw**  +l  "•  J*i"  *  ni  ^ 
p(a p)~  "  ®{  ***  +in*  Ai~A*  } 

+-|n*  ni.  ~  1 


P(m)-  aK-  '-§  »:  y  I nUl 
+B  -a^‘ + **  1 

F(s  4^  =  F  (-*3)  -  1 


F(i,0  =  Fl‘^)=  i 

-?&|.A»0‘  +  |y>.Jy.*-i/)‘{ 

Fl*,7)  -Ft7/*)  *  +  Ayl®  Ai~  7  J!‘  i 

f(4,4)=  a{/i?-aH 

F(4,5)=  F(s,4)  =  A- i  §  (*;-•«!  )} 

F  (4,4)  -?(<•,*)-  A-|  J  (aJ-/);  )} 
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F(4,7)  =  pU4)  -  A  {i  ( -1* -A*  )} 

F(5,s)  =  A {iW-At)]  * 

F(s>)= F(t,,s)  =  A  {f  (/>*-/)* ))  +G  Ll}i  ( nl-Ai )} 

fm =f(/,s)  -  a'  j  («;-<!)] + gl1  {f  ( *?-*? )} 

F(*>)«  AlilA,‘-A?)J  +  6L‘{|.(^;-Ai>} 

F(S7)  =  F^‘)*A}|«-flI>}  +  SL,-{ii(A.,-A,‘)) 

F(7,7)=  A{i(A;-7!*)}  +  6L1{|- 

C,  =  Cj  -  C3=  o 

-/***  (<-/£)  C5=-aj>ft^(/iJ-7i?)/3. 

c6  =  -j>cu(aJ-/£)/a.  c7=-*j><wU.s-<a- )  /s. 

A  computer  program  was  used  to  obtain  the  solution  of  the  problem.  A 
large  number  of  problems  were  solved  for  various  values  of  the  physical 
parameters,  r^ ,  r  ,  L  and  am,  and  the  propellant  properties  v  and 
E.  A  listing  of  the  program  Is  given  In  Appendix  A. 

DISCUSSION  OF  RESULTS 

It  has  been  noted  by  other  investigators  that  numerical  difficulties 
arise  in  the  analysis  of  propellants  when  Poisson's  ratio  approaches  0.5 
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special  variational  theorem  for  nearly  Incompressible  materials.  In 
trie  present  In/estlgatlon  the  numerical  problem  was  overcome  by  using 
a  CDC  6600  computer  with  double  precision.  This  results  In  approxi¬ 
mately  32  significant  digits,  sufficient  for  Poisson's  ratio  to 
0.4999999. 

A  study  was  conducted  on  the  affect  of  Poisson's  ratio  on  the 
propellant  stress  for  a  propellant  with  1.5  Inch  outer  radius,  0.5 
Inch  Inner  radius,  9.5  inch  length,  500,000  psl  Young's  modulus  and  a^ 
equal  to  8000  g.  Poisson's  ratio  was  varied  between  0.45  and  0.4999999. 
The  computer  program  results  for  or  at  the  propellant-Case  interface  at 
the  base  are  given  In  Table  1. 

Table  1 


V 

or(psi ) 

Av 

A  o 

e 

.45 

-1589 

-2.0 

X 

10‘3 

.48 

-1691 

.03 

-112 

CTi 

• 

00 

1 

X 

4 

10 

.49 

-1771 

.01 

-  80 

-2.1 

X 

_4 

10 

.499 

-2513 

.004 

-616 

-2.8 

X 

.5 

10 

.4999 

-3744 

•  .0009 

-1231 

-4.0 

X 

.6 

10 

.49999 

-4469 

.0009 

-725 

-4.8 

X 

_  7 

10 

.499999 

-4613 

.00009 

-144 

-5.0 

X 

.8 

10 

.4999999 

-4630 

.000009 

-  27 

-5.0 

X 

o 

1 

iO 
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MILLIMETER 


Figures  1  through  3  are  curves  showing  these  results;  the  first 
figure  being  a  direct  plot  of  stress  or  versus  Poisson's  ratio  v. 
Figures  2  and  3  show  more  clearly  that  the  stress  approaches  a  finite 
val«e  of  -4650  psi  as  Poisson's  ratio  approaches  0.5.  In  Figure  4  the 
dilatation  e  =  eQ  +  ef  +  ez,  at  r  «  rQ  and  z  =  0  is  plotted  against 
Poisson's  ratio.  The  resulting  curve  verifies  that  the  incompres¬ 
sibility  condition  e  approaches  zero  as  v  approaches  1/2  is  satisfied. 

In  addition  to  the  study  on  the  effect  of  Poisson's  ratio  on 
propellant  behavior,  analyses  were  made  to  determine  the  effect  of 
other  parameters  to  behavior.  The  results  obtained  show  that  the 
radial  stress  varies  linearly  with  respect  to  acceleration,  radii, 
length  and  Young's  modulus  of  elasticity. 

The  results  of  this  analysis  were  compared  with  two  other  sources; 
a  Rohm  and  Haas  finite  element  analysis  (ref.  3),  and  an  exact 
solution  for  an  incompressible  material  with  zero  Interior  radius. 

The  present  formulation  resulted  in  a  20  percent  higher  radial  stress 
than  the  finite  element  program,  and  2  percent  higher  than  the  exact 
solution,  at  the  propellant-case  Interface  at  the  base. 

CONCLUSIONS 

A  computer  program  was  developed  for  the  determination  of  the 
stress  distribution  due  to  a  slumping  propellant  contained  within  a 
launched  rigid  motor  case.  The  program  was  Initially  incorporated 
into  the  SATANS  finite  difference  program  for  the  calculation  of  the 
elastic  static  buckling  load  of  the  motor  case  subjected  to  high 
acceleration  loading.  The  stiffening  of  the  motor  case,  which  Is 
predicted  by  theory,  was  obtained  in  the  SATANS  analysis;  a  favorable 
structural  response. 
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An  unfavorable  effect  on  structural  Integrity  is  associated  with 
the  additional  contribution  to  yield  from  the  propellant  radial  stress. 
For  a  motor  case  of  2.5  Inch  outer  radius,  zero  Inner  radius,  1/8  Inch 
thickness,  20  Inch  length,  and  a  maximum  acceleration  of  8,000  g, 
the  radial  stress  of  approximately  10,000  psi  gives  a  hoop  stress  of 
100,000  psi.  The  contribution  to  yield  of  this  stress  component  is 
of  significant  magnitude.  That  is,  the  propellant  has  a  beneficial 
buckling  stiffening  effect  but  an  adverse  affect  on  yield. 

The  present  analysis  Is  based  on  the  linear  elastic  theory  of 
solids  and,  therefore,  the  results  are  not  valid  beyond  the  Initial 
yield  point  of  the  propellant.  Because  It  Is  difficult  to  obtain  the 
material  properties  of  propellants  (Young's  modulus,  Poisson's  ratio, 
and  Initial  yield),  definite  assertions  regarding  the  structural 
Integrity  of  the  prc reliant  are  not  given  here;  however,  the  stress 
results  obtained  In  some  of  the  analyses  show  that  additional  study 
regarding  this  concern  may  be  warranted. 
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STRESSES  IN  A  SLUMPING  FROPELLANT 
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